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We consider the X-ray properties of the redback class of eclipsing millisecond pulsars. These are transitional 
systems between accreting low-mass X-ray binaries and binary millisecond pulsars orbiting white dwarfs, and 
hence their companions are non-degenerate and nearly Roche-lobe filling. The X-ray luminosity seems to scale 
with the fraction of the pulsar sky subtended by the companion, suggesting the shock region is not much 
larger than the companion, which is supported by modeling of the orbital light curves. The typical X-ray 
photon spectral index is ~ 1 and the typical 0.3-8 keV X-ray efficiency, assuming a shock size on the order 
of the companion’s Roche lobe cross-section, is on the order of 10%. We present an overview of previous 
investigations, and present new observations of two redbacks, a Chandra observation of PSR J1628—3205 and 
a XMM-Newton observation of PSR J2129—0429. The latter shows a clearly double peaked orbital light curve 
with variation of the non-thermal flux by a factor of ~ 11, with peaks around orbital phases 0.6 and 0.9. We 
suggest the magnetic field of the companion plays a significant role in the X-ray emission from intrabinary 
shocks in redbacks. 


1. The Redback Population 


Millisecond pulsars are thought to be formed in bi¬ 
nary systems where an old neutron star is spun-up 
via long term accretion from an evolved companion. 
In recent years, the MSP recycling scenario has been 
dramatically confirmed through observations of so- 


called “redback” millisecond pulsar systems Roberts 


2011 which have non-degenerate companions and in 


some cases transition between states with no visible 
radio pulsations but with optical and X-ray evidence 
of an accretion disk, and a state where radio pulsa¬ 
tions are observed that regularly eclipse near supe¬ 
rior conjunction. The first of these transition objects, 
PSR J1023-I-0038, showed optical evidence for an ac¬ 
cretion disk in 2001 which had disappeared by 2004 
Thorstensen fc Armstrong|2005 . In 2007, radio pul- 


sations were discovered [Archibald et al.|2009 , and in 


2013 the M SP returned to an accreting state |Stappers 
et al.|[2M4 


Millisecond pulsars in compact binary systems have 
the potential of providing unique insights into pulsar 
winds. The companion forces a shock to occur at a 
distance ds only ^ 10^ times the light cylinder radius 
of the pulsar Ric = PscI2tt (where Ps is the spin pe¬ 
riod, and c the speed of light), as compared to the 
more typical ds ^ 10® — 10®R;c of the termination 
shock of pulsar wind nebulae around young, isolated 
pulsars. This means that the shock probes the wind 
in a region which might be significant in determining 
how the magnetization parameter a, the ratio of mag¬ 
netic energy to kinetic energy, goes from a presumably 


high value at the light cylinder to an apparently low 
value at the termination shock in typical pulsar wind 
nebulae [cf. Kennel fc Coroniti|1984 . The basic shock 
emission theory for such intrabinary shocks has gener¬ 
ally followed the outline of Arons & Tavani 1993 first 


developed for the original black widow system. In this 
model, the pulsar wind shocks with material ablated 
from the companion’s surface, which is presumably 
swept back around the companion and ejected from 
the system. In these models, it is generally assumed 
that the only significant source of magnetic field is the 
magnetization of the wind, and that the X-ray emis¬ 
sion is synchrotron which can be somewhat beamed 
either through a partially ordered magnetic field or 
doppler boosting. 

A Chandra observation of PSR J1023-I-0038 in its 
radio pulsar state revealed significant orbital variabil¬ 
ity over five consecutive orbits Bogdanov et al.|2011 


with a pronounced dip in the X-ray flux at supe¬ 
rior conjunction, when the companion is between the 
pulsar and observer and the intrabinary shock pro¬ 
duced through the interaction of stellar outflows is 
obscured. The X-ray spectrum consists of a domi¬ 
nant non-thermal component from the shock and at 
least one thermal component, likely originating from 
heated pulsar polar caps. The eclipse depth and du¬ 
ration imply that the shock is localized near or at 
the companion surface. However, the companion only 
subtends ~ 1% of the pulsar’s sky, so that if the wind 
is isotropic, only ~ 1% of the pulsar’s wind is inter¬ 
cepted by the companion, and only ~ 7% would be 
intercepted if the wind is confined to an equatorial 
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sheet. Bogdanov et al. 2011 inferred a high a from 
the estimated magnetic field of ^ 40G required to ac¬ 
count for the soft X-ray luminosity. 

An observation with NuSTAR of PSR J1023-I-0038 
just before it retu rned to the accreting state 


dulkar et al. 2014 showed that the spectrum of the 


intrabinary shock is a very hard power law (photon in¬ 
dex r = 1.17) with no apparent cutoff out to ~ 50keV, 
for a remarkable X-ray efficiency of ~ 2% of the total 
spin down power, or around all of the nominal spin 
down power in the wind that would be intercepted by 
the companion. Such a hard spectrum is not easily 
obtained from a pulsar wind nebula shock, and such 
efficiency is unprecedented. This might be an indi¬ 
cation of a significant equatorial enhancement in the 
wind or a significantly higher moment of inertia than 
the canonical lO'^^gmcm^, but it is still a remarkably 
high efficiency under any circumstances. 

Systematic studies of X-ray emission from redbacks 
show some commonalities. Linares 2014 examined 
the Swift XRT data on redbacks and noted that, while 
in the pulsar state, their 0.5-10 keV luminosities tend 
to be in the range of ~ lO^^ergs”^ divided into 
relatively high luminosity {Lx ^ lO^^ergs s“^) and 
relatively low {Lx ^ lO^^ergs s“^). Studies of indi¬ 
vidual systems show that, on average, there is orbital 
modulation with an overall increase of about a factor 
of 2 centered around inferior conjunction, often with 


a hint of a double peaked structure Bogdanov et al. 

2014|? 

Gentile et al. 

2014 

Kong et al. 

2012 . How- 


clearly distinguish fine structure to the orbital light 
curve. Black widows, on the other hand, show a much 
greater variety in their orbital light curves, with some. 


like the original black widow PSR B1957-I-20 Huang 


et al. 2012 , having peaks centered around superior 


conjunction and others around inferior conjunction 
Gentile et al. 2014 . On average, the redbacks are 


more luminous than the black widows in X-rays. 


The light curve modeling of Bogdanov et al. 2011 


suggests that the X-ray emission happens very close 
to the surface of the companion, which suggests that 
little of the wind that is not directly intercepted by 
the companion takes part in the X-ray emitting shock. 
The fraction of the pulsar’s sky subtended by the com¬ 
panion, He, can be calculated from knowledge of the 
relative masses (requiring knowledge of the orbital in¬ 
clination angle), the fraction of the Roche lobe the 
companion fills, and the orbital separation. The in¬ 
clination angle and Roche lobe filling fraction can be 
estimated from opt ical photometric light curves (eg. 
[Breton et al. 2013 ), and when combined with opti¬ 
cal radial velocity measurements and the pulsar orbit 
solution can be used to estimate the masses of the in¬ 
dividual components. In the table, we calculate Lie 
from our “best guess” estimates of neutron star mass, 
Roche lobe filling factor, and inclination angle using 
optical fits where available. On average, we estimate 
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Figure 1: Comparison of the X-Ray luminosity from the 
redback population to that of the black widow 
population. We define the shock luminosity as ELlc 
where Qc is the fraction of the pulsar sky subtended by 
the companion. We also plot the expected blackbody 
luminosity for each pulsar, assuming the relationship of 


Bognar et al. 


2015 
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Table I Redbacks 


Pulsar 

log .5“ 


r 

log Lx 


refs 

J1023-b0038 

34.7 

1.3 

l.UU_Q Qg 

32.0 

1.3% 

(1) 

J1227-4859 

35.0 

1.4 

-L--1-0_o.08 

31.9 

1.6% 

(2) 

J1628-3205 

34.2 

1.2 

1 0 + 0.8 
^‘^-0.7 

31.3 

1.1% 


J1723-2837 

34.7 

0.75 

1 1 0 + 0.02 
-*-•-*-^- 0.02 

32.1 

2.0% 

(3) 

J1816-b4510 

34.7 

4.5 

- 

31.0 

0.28% 

(4) (5) 

J2129-0429 

34.6 

0.9 


31.3 

1.2% 


J2215-b5135 

34.7 

3.0 

1 0+OA 
-‘-•^-0.3 

31.9 

1.4% 

(6) 

J2339-0533 

34.4 

0.4 

1 nQ+0-40 

l.Uy Qj^3 

30.6 

1.6% 

(7) ( 8 ) 


a. erg/s b. kpc, from dispersion measure except for 


J1023+0038 from parallax [Deller et al.|2012] and J1816+4510 
from optical [Kaplan et al.|p013| c. erg/s 0.3-8 keV d. 
estimated percentage of pulsar sky subtended by companion, 


(l)|Bogdanov et al.||2011| (2)|Bogdanov et aL| 

120141 (3) 

Bogdanov et al.||2014| (4)|Stovall et al.||2014| (5) 

Kaplan et | 

alJ 

20131 (6) 

Gentile et al.|2014| (7)|Romani Sz S 

o 

1 

(8) |Ray et al. 

|2014| 


Lie ~ 1.3% for redbacks and Lie ~ 0.3% for black wid¬ 
ows, accounting for much of the relative brightness 
of the shock emission of redbacks compared to black 
widows. 


We define a “shock luminosity” as ELle and plot 
that vs. the observed 0.3-8 keV X-ray luminosity 
of redbacks and black widows (Fi g0- We also plot 
the “expected” 0.3-8 keV blackbody emission from 
each pulsar based on a correlation determined from 
MSPs with good parallax meas urements log Lhh = 


(0.25±0.16)logF;-b(21.28±5.36) Bognar et al. 


2015 


We see that the shock luminosity and X-ray luminos¬ 
ity are correlated, with a typical soft X-ray efficiency 
relative to the shock luminosity of ~ 12%, albeit with 
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large scatter. We make no estimate of errors in the 
shock luminosity, being as they are dominated by the 
very uncertain distances in most cases and a lack of 
strong constraints from the optical data on inclination 
and the masses from the optical data in many cases. 
The redback with the smallest estimated flc and hence 
has one of the lowest luminosities is PSR J1816+4510. 
Optical studies of its companion suggest that it may 
be a proto-white dwarf which is significantly under¬ 
filling its Roche lobe Kaplan et al.||2013 


Spectrally, the X-ray emission tends to have a con¬ 
stant thermal component, presumably from heated 
polar caps and consistent with the typical thermal 
emission from MSPs, and an orbitally variable power- 
law component. The fit power-law tends to be very 
hard with photon spectral index P ^ 1 , harder than 
the typical spectra of pulsar wind nebulae around iso¬ 
lated young pulsars which have P ~ 1.5 in their inner, 
uncooled regions Kargaltsev fc Pavlov] 2010 . Below 
we report on new X-ray observations of two redbacks 
discovered by the Green Bank Telescope. 


2. PSR J1628-3205 

Discovered in a survey of Fermi sources with the 
GBT at 820 MHz (Sanpa-Arsa et al. in prep), PSR 
J1628—3205 is a 3.21 ms pulsar in a 5.0 hr or¬ 
bit around a companion with minimum mass > 
O.I 6 M 0 (assuming Mns = I. 4 M 0 ) (Hessels et al. in 
prep). The pulsar is eclipsed for about 20% of the 
orbit. It is modestly energetic with a standard spin- 
down energy of E = 1.8 x 10 ^'* ergs and an esti¬ 
mated distance from the pulse dispersion measure 
d ~ 1.2 kpc. Optical observations suggest it is Roche 
lobe filling with minimal heating of the companion [lI| 
et al.||2014 . 

We observed PSR J1628—3205 for 20 ks (slightly 
more than one orbit) on 05 May 2012 with the 
Chandra AGIS-S and detected ~ 180 counts. The 
counts as a function of orbital phase and energy, plot¬ 
ted in Figure suggest that there may be a dip 
in the above 2 keV flux near superior conjunction. 
The spectrum seems to have a significant power law 
component, with a purely blackbody spectrum not 
giving an acceptable fit. Using the CSTAT statis¬ 
tic of XSPEC (appropriate given the low number of 
counts per bin) suggests a pure power law fit pro¬ 
vides a somewhat reasonable fit (C-Statistic 24.97 
with 22 degrees of freedom), with best fit absorption 
nH = 1.3(0.2 — 2.5) x 10^^cm“^ and power-law index 
P = 1.60(1.23 — 2.00). Using the KS test statistic to 
determine goodness of fit results in 20 % of realizations 
having a lower test statistic, suggesting improvements 
can be made. Since most MSPs have a significant 
thermal component to their X-ray emission, we next 
tried an absorbed blackbody plus power-law fit. This 
resulted in a C-statistic of 21.08 with 20 degrees of 


J1 628—3205 Lightcurve 



J2215+5135 Spectral Fits 



Energy (keV) 


Figure 2: 20 ks Chandra ACIS-S observation of PSR 

J1628—3205. Top: Individual photon energies and 
average count rates as a function of orbital phase. The 
pulsar superior conjunction is defined as phase 1.25. The 
dashed lines show roughly the phase range of the radio 
eclipse. Bottom: 0.3-8keV spectrum showing absorbed 
power-law, blackbody, and blackbody-!- power-law fits. 


freedom, with less than 1% of KS realizations hav¬ 
ing a smaller test statistic. The best fit values were 
nH = 2.2 X lO^icm-^, kT = 0.20 keV and P = 1.14. 
The covariance between the blackbody temperature 
and the power-law index made it difficult to derive 
reasonable error bars if all parameters were allowed to 
vary freely, but by constraining the blackbody temper¬ 
ature to vary only between kT = 0.1—0.25 keV, within 
which range are the vast majority of MSPs, we find 
90% confidence regions of nH = (0.3—8.4) x 10^^cm“^ 
and P = (0.5 — 2.0). The 0.3-8 keV model flux is 
Fx = 8.8 X 10“^"^erg cm“^ s“^ with an unabsorbed flux 
of Fx = 1.2 X 10“^^ergcm“^ s“^, with roughly 70% in 
the power law and 30 % in the blackbody. T he fit nH 
is consistent with the [Drimmel et al. 2003 Galactic 
extinction model for a distance of 1.2 kpc. 
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Figure 3: XMM — Newton 0.1-10 keV light curve of 
PSR J2129—0429 as a function of orbital phase. The 
pulsar superior conjunction is indicated by vertical blue 
lines. 


PSR J2129-0429 



Figure 4: Unfolded XMM — Newton PN spectrum of 
PSR J2129—0429 at two orbital phases, keeping the 
blackbody component fixed. 


3. PSR J2129-0429 

Discovered in a survey of Fermi sources using 
the GET at 350 MHz Hessels et al.||2011|, PSR 


J2129—0429 is a 7.61 ms pulsar in a 15.2 hr orbit 
around a Me > O.37M0Companion which shows ex¬ 
tensive radio eclipses, as much as half the orbit at low 
frequencies (Hessels et al. in prep). The pulsar has a 
very high magnetic field for a MSP (H ~ 1.6 x 10® G), 
and so still has a high spin down energy E ~ 3.9 x 10®"^ 
despite its relatively long spin period. The dispersion 
measure distance is d ~ 0.9 kpc. A variable, bright 
UV counterpart was evident in the Swift UVOT, as 
was significant X-ray variability from the Swift XRT 
data. Further optical observations suggest the com¬ 
panion is minimally heated and mostly Roche lobe 
filling and radial velocity measurements suggest a 
pulsar mass > 1.7MQand a companion mass 


Me ~ 0.5Mq Bellm et al. 2013 . These system prop¬ 


erties suggest that PSRJ2129—0429 is in a relatively 
early stage in its evolution compared to other redbacks 
which are more fully spun-up and have typical mag¬ 
netic fields of a few 10® G. Very large orbital variations 
are observed through radio timing, and pulsations are 
dominant in the y-ray emission. 

We observed PSR J2129—0429 for 70 ks with 
XMM — Newton. There were no background flares 
during the observation, meaning we got continuous 
coverage over slightly more than a complete orbit. 
The X-ray light curve has very large amplitude vari¬ 
ations, with two clear peaks centered on the pulsar’s 
inferior conjunction (Figj^. We first fit the spectrum 
with an absorbed blackbody plus power-law, which 
gave an adequate fit. The flux is dominated by the 
power-law component, with an average 0.3-8 keV flux 
Fj: = 2.25 ± 0.05ergcm“® s“^. There is very little 
absorption {nH = 1.8(0 — 4.6) x 10®®cm“®) and the 
thermal component {kT = 0.21(0.16—0.26) keV) has a 
0.3-8 keV flux F^b 1.2x lO^^'^ergcm”® s“^, or about 


1/4 of the flux near superior conjunction. The power- 
law component is very hard (F = 1.04(0.92 — 1.15)), 
similar to other redbacks. Presuming a constant ther¬ 
mal component throughout the orbit, the difference 
in the non-thermal flux between the peak at orbital 
phases 0.575-0.65 and the minimum at phases 0.2-0.3 
is about a factor of 11 (Fig. |^. There is no evidence 
of significantly increased absorption. Complete spec¬ 
tral results will be presented in an upcoming paper 
(Roberts et al. in prep). 

This very remarkable variability suggests that a 
large fraction of the shock region is blocked by the 
companion around superior conjunction, suggesting a 
quite small emission region and a large inclination 
angle. The two distinct peaks may be a result of 
doppler boosting and/or relativistic beaming of the 
synchrotron radiation. The latter would require a 
strong, well ordered magnetic field. The orbital phases 
of the peaks, ~ 0.6 and ~ 0.9, are quite curious. If 
the shock was wrapped around the companion, then 
you would expect there to be peaks between phases 
0.0-0.5. The qualities may suggest a significant role 
for the magnetic field of the companion. If the com¬ 
panion is tidally locked, like one would expect, then 
the orbital period of 15.2 hr is the spin period of the 
companion, which is very rapid. Low mass, rapidly 
spinning stars can have surface magnetic fields of sev¬ 


eral hundred to a few thousand Gauss Morin 2012 


Such potentially large companion fields should not be 
ignored when investigating the shock emission from 
redbacks. 

In summary, X-ray emission from the intrabinary 
shock in redbacks is orbitally dependent, with the in¬ 
creased emission centered on inferior conjunction with 
potentially a fairly ubiquitous double peaked struc¬ 
ture. The emission seems to come from a region that 
is not much larger than the companion, is very hard 
and very efficient, which needs explanation. The pre¬ 
viously ignored potential role of the companion’s mag- 
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netic field in the shock dynamics needs to be consid¬ 
ered. 
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